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Scanning tunneling spectroscopy (STS) enables the local, energy- resolved investigation of a samples surface 
density of states (DOS) by measuring the differential conductance (dl/dV) being approximately 
proportional to the DOS. It is popular to examine the electronic structure of elementary samples by 
acquiring dl/dV maps under constant current conditions. Here we demonstrate the intricacy of STS 
mapping of samples exhibiting a strong corrugation originating from electronic density and local work 
function changes. The confinement of the Ag(l 1 1) surface state by a porous organic network is studied with 
maps obtained under constant-current (CC) as well as open-feedback-loop (OFL) conditions. We show how 
the CC maps deviate markedly from the physically more meaningful OFL maps. By applying a 
renormalization procedure to the OFL data we can mimic the spurious effects of the CC mode and thereby 
rationalize the physical effects evoking the artefacts in the CC maps. 

Shortly after the initial demonstration of scanning tunneling microscopy (STM) 1 , theoretical analysis of 
electron tunneling through the vacuum barrier has shown that the tunneling current I T depends on the local 
density of states (LDOS) of the sample 2 ' 3 . Subsequently, scanning tunneling spectroscopy (STS) was intro- 
duced, providing access to the energy dependence of the LDOS 4,5 and enabling the local characterization of 
physical and chemical properties of conducting samples. The powerful capabilities of STM/STS have been 
demonstrated in the investigation of the local electronic properties of surfaces 6 " 11 . The key advantage of STS 
over the classical, spatially averaging spectroscopy techniques is its capability to investigate the electronic struc- 
ture on the atomic scale, while at the same time providing topographic information. Therefore, STS is perfectly 
suited to study adsorbates positioned on substrates 12 " 15 and confinement effects of surface state electrons 16 " 24 , 
providing a complementary approach to /c-vector resolving photoemission spectroscopy techniques 24 " 26 . 

STS is typically carried out by measuring the differential conductance (dl/dV) under open feedback-loop 
(OFL) conditions by applying to the bias voltage V B a small harmonic modulation with a specific frequency 
oj mo d and extracting the corresponding modulation -related signal via a lock- in technique. Within certain approx- 
imations 27,28 , the recorded dl/dV signal is proportional to the LDOS p s (r 9 E) of the sample at the position of the 
STM tip r and energy E. When focusing on the spatial distribution of the LDOS of a specific area (x, y), 
spectroscopy data is often acquired by recording dl/dV while simultaneously scanning the tip with a constant 
I T in order to obtain contour maps for a certain energy E B = eV B . However, as was already demonstrated ear- 
lier 29 " 36 , in this so-called constant-current (CC) mode, the dl/dV signal does not only depend on the DOS, but is 
convoluted with the change in tip -sample separation z. This convolution results from the dependence of I T and 
dl/dV on the transmission probability T as can be seen when expressing the quantities within the 1D-WKB 
approximation 7,27 ' 36 ' 37 . Since T depends exponentially on z both quantities change drastically upon variation of z. 
In simple systems meaningful LDOS information can be retrieved from the convoluted CC data by estimating the 
change of T from the z-signal 36,38 . However, when complex systems with appreciable surface corrugation and local 
variation of the work function are under investigation, T depends non-trivially on z and this method fails. 

In this work we first shortly assess the theoretical background for space-resolved dl/dV measurements of 
inhomogeneous samples. Then, we discuss the effect of a varying tip-sample separation during a two-dimensional 
(2D) measurement on the experimental spectroscopy data of molecular networks confining the quasi-free 2D 
electron gas of the Ag(lll) surface state to periodic quantum dot arrays 22,23 . 
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By comparing datasets of 2D dl/dV maps acquired in both CC and 
OFL mode, we directly relate the different measurement techniques 
and exemplify the strong deviations of the CC maps from the phys- 
ically straightforward OFL counterparts. We thus demonstrate that 
the CC data can be misleading even on a qualitative level for the 
investigated corrugated systems. Finally, we introduce a special 
renormalization procedure, accounting for the locally changing T 
which allows us to reproduce the CC maps from the OFL dataset. 
The successful reproduction is consistent with the theoretic under- 
standing of the spectroscopy modes. 

Results 

Theoretical considerations. In a first step, we review the dependence 
of the STS signal on the transmission probability and the density of 
states Ps(E) between the Fermi energy E F = 0 of the sample and the 
energy E B = eV B at which the dl/dV signal is measured. In the WKB 
approximation for planar electrodes 6 I T depends on the DOS of the 
tip, p T , the DOS of the sample, p s , and on the transmission 
probability T: 

eV B 

-E)p s (r T ,E)T(z,eV,E)dE (1) 
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In this description, the transmission probability for electrons with 
energy E to tunnel across the vacuum barrier is given by 
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with the tip-sample separation z(r) at the lateral position r = (x,y) 9 
the effective electron mass m, and the work functions of the tip and 
the sample materials, <X> r and O s respectively. 

Even when p T is assumed to be constant, the dl/dV signal is no 
longer directly proportional to p s when T changes perpetually while 
dl/dV is measured in CC mode: 

eV B 
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In this expression, the second term accounts for high energy con- 
tributions which we neglect in this consideration. Instead, we con- 
centrate on the influence of p s and T to account directly for the 
experimental conditions when acquiring dl/dV maps in CC mode 
at low bias. For calculations we therefore use the expression 

d/ azp s (eV B )T(z,eV B ) (4) 

V B 



dv 



for the differential conductance. 

Origin of artefacts in the spectroscopic data recorded in the CC 
mode. To illustrate the difference between dl/dV measurements in 
OFL and CC mode, we first discuss three exemplary cases of E- and x- 
dependencies (gray, blue, and red in Fig. 1) in a simple model system, 
namely a topographically featureless, flat surface, divided into two 
regions I and II being traversed during the measurement. Assuming a 
W tip over the Ag(lll) surface, the work functions were chosen as 
<I> r = 4.55 eV and O s = 4.74 eV and m = 0.42 m e was set for the 
Ag(lll) surface state electrons 25 ' 39 . Without loss of generality we 
focus on the energy range from E F = 0 eV to E x = 0.1 eV. 

In the first example (gray with squares in Fig. 1), we consider two 
different energy dependencies of p s (E, x) in the two regions I and II 
(Figs, la and b) resulting in the same state density at the energy of 
interest £ l5 namely p s (E, x<l) increases from 0 for E = 0 to 1.5 for E 
= Ei and Ps(E, x > 1) decreases from 3 to 1.5 in the same range. Thus, 
at E = Ei the state density is the same in both regions and accordingly 
the measurement of ps(Ei) along x should give a region independent 
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Figure 1 | Model system and simulated dJ/dVtraces. (a) and (b) Different 
energy dependencies for the local density of states ps(E) are assumed in 
regions / and II. (c) For the two regions different spatial distributions for 
p s (xy Ei) at energy E x represent the reference quantities that are to be 
reproduced, (d) Simulated dl/d V traces along x obtained in CC mode. The 
simulated traces deviate strongly from the reference signals in (c). In the 
red (upward-pointing triangles) and gray (squares) scenario, the p s (x) 
features are not reproduced correctly, in the blue case (downward- 
pointing triangles) the step is reproduced only qualitatively. 



constant signal (Fig. lc). However, the calculated CC mode dl/dV 
signal (Fig. Id) exhibits a step. The CC signal was simulated by first 
finding z and Tby self- consistently solving Eq. (1) and (2) for ident- 
ical I T in both regions and then calculating the dl/dV signals in both 
regions via Eq. (4). The CC mode dl/dV signal is smaller in region II 
because the integrated density of states between E = 0 and E = E x is 
larger in region II causing a retraction of the tip for establishing the 
constant I T and thereby a reduction of T. 

In a second and third example (red with upward-pointing trian- 
gles and blue with downward-pointing triangles in Fig. 1) we con- 
sider an identical step- wise increase in p s (x, E{) when moving the tip 
from region I to region II (Fig. lc), combined with different linear 
energy dependencies of p s (x, E) in the two regions (Fig. la and b, red 
and blue). The calculated dl/dV signal (Fig. Id) exhibits an exagger- 
ated step (factor 2) in the blue case and an inverted step in the red 
case. 

These examples demonstrate the strong impact of T on the dl/dV 
signal obtained in the CC mode where the tip-sample separation is 
perpetually adjusted to maintain a constant I T . Thus, the spatial 
dependence of a dl/dV signal recorded at an energy E x in CC mode 
can only in special cases be associated with the actual spatial depend- 
ence of the sample LDOS at that energy. In general, such mea- 
surements suffer from incorrect intensity ratios and yield even 
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qualitatively wrong behavior when comparing spectral densities at 
different positions. 

The artifacts in the dl/dV signal resulting from a spatially non- 
constant T can not only result from changes in the tip -sample gap 
distance z as discussed so far, but also from local variations of the 
sample work function <I> S , even when the tip is scanned at a constant 
height. Therefore, we address the effects of local O s variations on dll 
dV measurements in CC mode. To this end we assume a flat, fea- 
tureless surface with constant p s (x, E) and with two regions, one 
having a low work function c|) Sjl = 3.74 eV and one having a <E S;2 
= 4.74 eV. When we calculate the CC-dJ/dy signals, the deviation in 
both regions is quite low. Even though a reduced work function 
results in a higher tip-sample separation, the integrated transmission 

reV 

probability TdE is identical (cf. Fig. 2). However, the difference 
Jo 

for the absolute values of T at the setpoint (here: I T = 0.1 nA, E B = 
0.2 eV) is notable (1.5%). Thus, care has to be taken, when dl/dV 
spectra are measured on a system with variable <I> S . 

Imaging a porous molecular network. In order to obtain a 
physically meaningful comparison of spectral intensities at 
different positions on the sample, we use the following procedure 
which should in principle be applicable to any other system. It is 
important to carefully select a proper setpoint for stabilizing the tip 
before recording the OFL spectra. In our case of the confined surface 
state on Ag(lll) we choose E B = -0.2 eV and I T = 0.1 nA. At this 
energy the contributions of the molecular orbitals are small and the 
surface state-related DOS is not present, i.e., the sample DOS 
originates from the projected bulk DOS of Ag only. Accordingly, 
the dl/dV signal is to a good approximation constant in the whole 
area (x,y) of interest. It follows that the artifacts produced by local 
tunneling gap changes originating from local work function changes 
can be eliminated by normalizing all spectra such that they have the 
same dl/dV value at the setpoint energy, see Eq. (4). This is, because 
under such tunneling conditions, the tip in CC mode follows a 

reV 

contour of constant Td£, which can be compensated for, in a 
Jo 

first approximation, by the mentioned normalization condition. In 
our case, where the dl/dV is featureless in the energy range between 
—0.2 V and —0.1 V, the arithmetic mean value m* of all the bias 
points within this energy interval of the dl/dV was calculated for each 
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Figure 2 | Calculated T in two otherwise featureless sample regions with 
different work functions <I>. Whereas the integrated transmission 
probability which contributes to I T is identical in both cases (low O 
displayed in red with upward-pointing triangles, high O in blue with 
downward-pointing triangles), the absolute values differ notably (1.5%) at 
the selected setpoint (I T = 0.1 nA, E B = 0.2 eV). 



spectrum i. Subsequently, each spectrum was multiplied by a 
normalization factor = f/m b where / is in principle a deliberate 
value, but typically the mean value of all m z is used. This allows us to 
directly compare spectral intensities within the pores and on the 
molecules in a unified framework. A further normalization via the 
division by the measured II V, which is typically performed for 
semiconducting surfaces 7,40 , was not performed due to the 
approximately linear II V encountered on the metal surface in the 
energy interval of interest. 

In order to understand and qualitatively assess the challenges 
addressed above, we study the effects of CC mode dlldV measure- 
ments in the model system of Ag(l 1 1) 2D surface state electrons that 
are confined to 0D quantum wells via organic molecular networks 
self-assembled from sexiphenyl dicarbonitrile (for structure model, 
see Fig. 3) 41 . The confining properties and resulting bound states 
were described experimentally and assessed via Boundary Element 
Method simulations in previous studies 22 ' 23 . 

The exact contrast mechanisms for molecules being imaged in 
STM is currently still widely discussed 42 " 47 . Molecules often markedly 




Distance along line [nm] 

Figure 3 | Imaging parameter dependent appearance of an organic 
network. Top: structure model of sexiphenyl dicarbonitrile (C black, 
N blue, H white), (a) One possible network among the different self- 
assembled, organic phases is the kagome lattice. The open-porous network 
follows a four-fold bonding motif and each molecule interacts with five 
adjacent molecules (inset). I T = 0.1 nA, V B = —0.5 V (b) Height profile 
extracted along the red line in (a) shows that the apparent molecule height 
at these tunneling parameters is ~ 1.2 A (red profile). Depending on the tip 
and tunneling parameters, the molecule height varied between 0.2 A and 
1.6 A (black dashed lines display upper and lower limit). The blue profile 
shows the calculated tip-sample separation when a local lowering of the 
electron work function by 1 eV at the position of the molecules is assumed 
at I T = 0.1 nA and V B = 0.2 V. 
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appear in STM topographs at bias voltages in the gap between the 
lowest unoccupied and highest occupied molecular orbitals where 
the orbital-related contributions to the tunnel current should be 
irrelevant. Therefore different effects are necessary to explain this 
phenomenon. Among others, a local variation of the work function is 
a possible candidate to rationalize molecule imaging. In the case of 
self- assembled, organic sexiphenyl dicarbonitrile networks in which 
the molecules appear as stick- like protrusions and the network bond- 
ing follows a four-fold nodal motif, one of the possible geometries is 
the kagome network (see Fig. 3a). When a topographical line-scan 
across the molecule network is extracted along the red line in Fig. 3a, 
the apparent height of the molecules varies depending on the experi- 
mental conditions in the interval of —0.2-1.6 A. In Fig. 3b the actual 
line scan (red, obtained at V B = —0.5 V) is displayed along with the 
upper and lower limit of the observed height (dashed, black lines). 
From STS we know that the LUMO of these molecules is positioned 
at around V B = 1.5 V. Taking into account that for polyphenylene 
molecules the HOMO -LUMO gap is typically more than 3 eV, the 



line scan was recorded near to the center of the HOMO -LUMO gap 
suggesting that the contribution of the frontier orbitals to the 
molecular contrast is very weak. The deposition of sexiphenyl mole- 
cules onto the Au(lll) surface at monolayer coverage results in a 
lowering of the electron work function by 0.73 eV 48 . In the model 
system presented above a local reduction of <I> S by 1 eV at constant p s 
results in the tip being retracted by —0.5 A. This is in in the same 
ballpark as the apparent molecule height of —0.2-1.6 A in CC topo- 
graphs and supports our reasoning that a significant molecular cor- 
rugation can be explained solely by work function considerations. 

Comparing spectroscopic information obtained in OFL mode vs 
CC mode. When Ag(l 1 1) surface state electrons are confined within 
the molecular network, a set of localized electron states evolves. We 
first probe the system via OFL point spectroscopy (see Fig. 4). When 
we investigate specific positions within the network pores (Fig. 4a), 
we can identify the corresponding LDOS signatures (Fig. 4b, spectra 
in the quasi-hexagonal pore are offset from the spectra in the 
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Figure 4 | OFL dJ/d V spectra for a self-assembled supramolecular kagome network on Ag(l 1 1). (a) In the topograph the positions of OFL point spectra 
are indicated. Colored symbols mark individual spectra displayed in (b) and the red line defines the positions of a total of 50 spectra visualized in (c). 

(b) At selected positions in the pores of the networks, different dl/dV fingerprints are identified (spectra in the two pore types are offset for clarity). 

(c) Color plot of 50 OFL spectra taken along the line indicated in (a) revealing a strong spatial variation of the LDOS within the network cavities. The 
expression of quantum confinement leading to defined electronic states becomes clear. The selected positions reproduced in (a and b) are indicated by 
respectively colored dashed lines, (d) A dl/dVprofile (red solid line) extracted at —50 mV from the set of OFL spectra (displayed in (c) with the white 
dashed line) clearly shows the onset of the surface state in both pore types. On top of the molecules (represented by the peaks of the topographic line drawn 
with the black dashed line) the surface state is suppressed, i.e., low dl/dV intensity is found in agreement with their scattering properties, (e) In CC mode a 
higher, nearly constant spectral signal (red solid line) is observed across the smaller, triangular pores and molecules, whereas the quasi-hexagonal pore in 
the center shows a minimum. Upon renormalization in the energy interval [-50 mV. . .0 mV] the OFL data (blue solid line) reproduces the observations 
of the CC mode (related topography profile indicated by a black dashed line). 
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triangular pores for clarity). In the center of the quasi-hexagonal pore 
(blue) two features are observed at —5 mV and 280 mV, 
respectively. Towards the cavity rim (red and green) the intensity 
of the peaks at —5 mV and 280 mV decreases, however an 
intermediate feature at 80 mV is found. In the smaller triangular 
pores (black and magenta) the first state is shifted towards higher 
energies (70 mV) due to the reduced quantum dot size. 

When all 50 spectra acquired along a line across the different pores 
(Fig. 4a, red line) are depicted in a color plot (Fig. 4c, respectively 
colored dashed lines indicate positions of point spectra displayed in 
Fig. 4a), the spatial and energetic variation of the LDOS becomes even 
more clear. In the central cavity the states featuring 0, 1 and 2 nodes 
can be identified. The triangular pores show an energetically broader 
intensity maximum, with the expression of individual states being less 
clear. A deeper analysis of the system was presented in Ref. 22. 

The kagome model system will now serve for a comparison of the 
CC and OFL spectroscopy techniques. In a first step, we compare dll 
dV signals acquired across the open pores of the self-assembled 
network, see Fig. 4a. At E B = —50 mV the dl/dV intensity extracted 
from the 50 OFL point spectra presented above (Fig. 4d) is related to 
measurements performed in CC mode (Fig. 4e, red). Distinctively 
different dl/dV distributions are measured in both operation modes. 

The OFL technique leads to maxima at the center of the pores, 
whereas on top of the molecules minima are found (Fig. 4d, red). In 
contrast, dl/dV measurements at identical energy in CC mode 
contradict this finding (Fig. 4e, red). Here, the central pore shows 
a depression, whereas the molecules and smaller cavities show 
approximately constant dl/dV at higher intensity. Thus the 
LDOS variation obtained in CC mode fails to qualitatively repro- 
duce the behavior expected for confined electrons in the network 
cavities in this strongly corrugated system, even relatively close to 
E F . However, the influence of T on the CC mode can be reproduced 
by employing the normalization procedure presented earlier in the 
energy interval from E F to E B to the OFL dataset instead of the 
spectroscopically featureless setpoint region, yielding the blue curve 
in Fig. 4e. Via this renormalization the dl/dV signature observed in 
CC mode is qualitatively reproduced by the renormalized OFL 
dataset. 

When we extend the comparison of CC and OFL dl/dV measure- 
ments towards 2D mapping of the LDOS, the strikingly different 
observations for both modes in the line- scan are further elucidated. 
We directly relate maps acquired in OFL and CC mode at — 50 mV 
and 250 mV in Fig. 5. The OFL dl/dV maps were extracted from 
point spectra on a 84 X 84 pixel grid equally distributed over the area 
displayed in the topographic image in Fig. 5a, where for each spec- 
trum the tip was stabilized at E B = —200 mV and I T = 0.1 nA. 
Figures 5b and e represent the OFL maps for E B = —50 mV and 
250 mV respectively and are given with the identical color code as in 
Fig. 4c. In comparison, CC dl/dV maps are shown in Fig. 5c and f. 

The —50 mV OFL map (Fig. 5b) is in accordance with the line of 
spectra shown in Fig. 4. As expected, we observe the onset of the first 
bound state within the triangular and the quasi-hexagonal pores, 
whereby the molecules show the lowest intensity. The CC dl/dV 
map at the same bias voltage (Fig. 5c) shows increased spectral den- 
sity above the molecular backbones and triangular pores, whereas the 
quasi- hexagonal pores show low intensity. Note that the CC map was 
recorded over a domain exhibiting the opposite chiral enantiomorph 
of the network which is barely visible in the features of the network, 
but explains the different positions of the centers of the pores. From 
our earlier work 22 we know that the chirality signature does not 
influence the other properties of the confined electron states. When 
the OFL dataset is renormalized in an identical fashion as in Fig. 4e, 
the dl/dV map changes drastically. The renormalized OFL map 
(Fig. 5d) now closely resembles the data acquired in CC mode. The 
maximum intensity is also found on the molecules and triangular 
pores, whereas the quasi-hexagonal pores show depressions. Thus, 
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Figure 5 | Comparison of 2D dJ/d Vmaps acquired in CC and OFL mode. 

(a) Topograph image of the area over which the OFL spectra were 
distributed (I T = 0.1 nA and V B = -0.2 V). (b)-(d) E B = -50 mV: 

(b) The dl/dV map obtained in OFL mode displays the onset of the first 
confined state in the triangular and quasi-hexagonal pore, respectively. 

(c) The map acquired in CC mode shows higher intensity over the 
molecular ligands and triangular pores, whereas the quasi-hexagonal pores 
feature depressions, (d) After renormalization of the OFL map over the 
energy range E F ^ E B , the observations found in CC mode are reproduced 
qualitatively, (e)-(g) E B = 250 mV: (e) The OFL map shows confined 
states in both pores alongside depressions at the positions of the molecules, 
(f) The CC map exhibits high intensity on the binding motifs and in the 
center of the quasi-hexagonal pores. Thus the spatial dependence of Ps(E B ) 
is obtained in an even qualitatively wrong manner, (g) The renormalized 
OFL map reproduces the CC map in great detail. Relative intensities and 
spatial distributions of the observed features show good agreement. The 
color-code of all maps is identical to that used for Fig. 4c. 

a good agreement between the CC dl/dV map and the renormalized 
OFL map is achieved. 

Further away from E F , at E B = 250 mV, the OFL dl/dV map 
shows different confined states in both pores (Fig. 5e). In the 
quasi-hexagonal pore, maxima are found in the center and close to 
the corners. The triangular pores exhibit a density distribution that 
resembles a three winged propeller. In comparison, the CC dl/dV 
map again differs markedly from the expected distribution of con- 
fined electron states (Fig. 5f). The maximum intensity is now found 
on the binding motifs of the molecular network, leading to six rect- 
angular features surrounding the quasi-hexagonal pore. In the center 
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Figure 6 | Effect of partial renormalization intervals, (a) OFL map 

acquired at E B = 300 mV. (b) CC map acquired at E B = 300 mV. When the 
OFL map is renormalized only over either half of the entire energy interval 
(E F ^> 150 mVin (c) and 150 mV^ 300 mVin (d)) the relative intensities 
of the different features and and their related spatial distribution are only 
poorly reproduced, (e) When the OFL data are renormalized in the energy 
interval E F ^> E B the obtained map agrees quantitatively with (b). 

we observe a maximum, surrounded by a hexagonal depression. In 
contrast, the triangular pores feature minimal intensity. Therefore, in 
CC mode merely the maximum in the center of the large pore is 
reproduced. After the OFL map was renormalized with the normal- 
ization interval ranging from E F to 250 mV, the map nicely matches 
that acquired in CC mode. The overall spatial distribution as well as 
the relative intensity of the observed features agree well. 

Discussion 

The analysis presented exemplifies that by renormalizing the OFL dll 
d ymaps, the data from CC dl/dV maps can be reproduced in great 
detail (see the Supporting Information for an extended series of dll 
d ymaps). This effect can be attributed to the operation of the feed- 
back loop in CC mode. When I T is kept constant for a given bias 
voltage, the tip-sample separation z is perpetually adjusted to com- 
pensate for changes in p s during the scan. Therefore, by maintaining 
a constant I T , the feedback loop operation entails an integration over 
the LDOS between E F and E B to account for the spatial changes in 
tunneling probability. By employing the renormalization procedure, 
this effect is reproduced mathematically from the OFL spectra. When 
all spectra are subsequently renormalized to their mean dl/d lvalue 
in the interval E F E B > the LDOS contribution to the tunneling 
current and thereby the readjustment of T at every position is simu- 
lated. The relative intensities of all spectra are shifted depending on 
the ratio between the mean spectral intensity and the actual intensity. 
Thereby, the effect of a variable tip-sample separation and finally the 
dl/dV maps acquired in CC mode are reproduced. 

Finally we demonstrate that all states in the range between E F and 
E B contribute to the artifacts obtained in CC mode. In Fig. 6b, a CC 
dl/dV map acquired at E B = 300 mV is shown alongside the renor- 
malized (E F — > E B ) OFL map in Fig. 6e. At this energy, the rectangular 
features located on the binding motifs and the maximum in the center 
of the quasi-hexagonal pore exhibit an identical intensity in the CC 
dl/dV map and the renormalized OFL map. When the OFL map is 
renormalized over the energy ranges E F —> 150 mV and 150 mV —> 
E B (see Fig. 6c and d), respectively, the maps are already similar to the 
CC map. However, the relative intensities are qualitatively wrong. In 



Fig. 6c the rectangular units show higher intensity than the maximum 
in the center of the quasi-hexagonal pore. Furthermore, the peak in 
the center has a smaller diameter than in the CC dl/dV map. 
Contrarily, in Fig. 6d the central maximum is spatially more extended 
with higher intensity than the rectangular units on the binding motifs. 
Thus, only the renormalization over the entire energy range can 
reproduce the findings from CC dl/dV mapping in detail. 

In summary, we conducted a detailed STS investigation of the 
surface state confinement in a nano-porous organic network focus- 
ing on the comparison of spectroscopic data recorded in both con- 
stant-current and open-feedback-loop mode of the very same 
system. The spurious effects of the variable tip -sample separation 
during CC measurements on highly corrugated samples with lateral 
work function variations were reported and analyzed. Our findings 
demonstrate and explain pitfalls that frequently occur in this method 
and thereby elucidate the intricacy of mapping the LDOS in two 
dimensions with STS. Especially for inhomogeneous systems great 
care has to be taken experimentally and during data processing. 
While CC dl/dV mapping is a viable tool to explore spatial distri- 
bution of electronic states in simple samples, physically meaningful 
LDOS maps of inhomogeneous systems can only be obtained when 
OFL spectroscopy is employed. Albeit acquiring datasets of OFL 
spectra on a grid is time consuming and experimentally challenging, 
it delivers a complete analysis of the energy range under investigation 
as opposed to recording individual CC maps. On the basis of a well- 
defined setpoint, OFL datasets enable a physically meaningful com- 
parison of spectral intensities at different energies and positions. It is 
suggested that such considerations are generally relevant for the 
detailed understanding of electronic properties in organic-inorganic 
hybrid systems. 

Methods 

Scanning tunneling microscopy/spectroscopy. STM measurements were 
performed using a home-built LT-STM 49 50 in UHV (3* 10" 11 mbar base pressure) at 
temperatures below 10 K. The Ag(lll) surface was prepared by repeated cycles of 

Ar + sputtering (flux— 7-^— -) and annealing to 740 K, in order to obtain atomically 
cm 2 

flat terraces separated by monatomic steps. Sexiphenyl dicarbonitrile molecules 
(for structure model, see Fig. 3) were deposited from a quartz crucible in an organic 
molecular beam epitaxy (OMBE) source at 585 K, while the substrate was held at 
room temperature. All topography was acquired at constant tunneling current with 
closed feedback loop. dl/dV signals were recorded with bias modulation at ~ 1400 Hz 
at a modulation amplitude of 5 mV rms, while the lock-in time constant was set to 
either 20 or 50 ms. 
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